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Abstract 



The problem of biped locomotion at steady speeds is discussed through 
a Lagrangian formulation developed for velocity-dependent, body driving 
forces. Human walking on a level surface is analyzed in terms of the data on 
the resultant ground-reaction force and the external work. It is shown that 
the trajectory of the center of mass is due to a superposition of its rectilin- 
ear motion with a given speed and a backward rotation along a shortened 
hypocycloid. A stiff-to-compliant crossover between walking gaits is de- 
scribed and the maximum speed for human walking, given by an instability 
of the trajectory, is predicted. 
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In biology, many fundamental discoveries come from studies of animal movement. The 
integrative approach to the locomotion of animals focuses on the interaction between the 
muscular, tendon and skeletal subsystems and the environment. [1] 

During locomotion, the body as a whole performs several functions. Chemical energy 
released by muscles and mechanical elastic energy stored in passive muscles and tendons [2] 
are transformed into external and internal work, [3] and are partially lost as a heat. The 
muscular system provokes the ground-reaction forces applied to the animal body. The resul- 
tant force, including gravity and air resistance, accelerates and decelerates the body's center 
of mass (COM). This body driving force is therefore involved in level walking and running, 
even when the average velocity remains constant. Over a complete step cycle, the driving 
force performs a certain external work to maintain a given speed. Studies of the mechanical 
efficiency [3] of animal locomotion at different steady speeds, determined through the external 
work and oxygen consumption, provide evidence that walking is more energetically econom- 
ical than running. This finding corroborates the old idea that walking in humans, primates, 
and ground-dwelling birds can be understood as swings of an ideal pendulum. Indeed, the 
body vaults up and over each leg in an arc in each step, similarly to an inverted plane pendu- 
lum, and kinetic energy is transformed into gravitational energy when the body falls forward 
and downward. Such a stiff-legged mechanics of walking modeled by the compass-arc in- 
verted pendulum is widely employed, [1,3,4] but by no means exhaustive. [1,5] Unlike the 
swing pendulum dynamics driven by constant gravitational force, the dynamics of animal 
biped walking (and quadruped trotting) is accompanied by body undulations, pulses, and 
peristaltic waves [1] generated by the ground-reaction force. In this Letter, we employ the 
fundamental principles of classical mechanics to approach to the problem of level locomo- 
tion. At very low speeds, animal movement is treated through the linear vibrations of the 
body near its quasistatic equilibrium given by quiet standing. [6] Nonlinear body motion 
effects, controlled by velocity-dependent forces, are taken into account at higher speeds. 

All the three components (forward, vertical and lateral) of the resultant force applied 
to the ground are measured with good accuracy by means of the force-platform techniques. 
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[7] The lateral body displacements are relatively small and, thus, the COM motion can be 
fairly described by the instantaneous polar vector R(i) defined in the ground coordinate 
system (see Fig. 1). For walking at a steady speed V, it is convenient to exclude the 
translational degree of freedom by introducing r(t), with x(t) = Vt and y(t) = H, where H 
is the height of the COM. In that way, the libration motion is given by Ar(t) = R(t) — r(t). 
The corresponding driving force AF(t) follows from the force-platform records: the ground- 
reaction force F(t) is observed as the oscillating force near the body weight, thus, F(t) = 
-mg + AF(t). Taking into account that the muscle-tendon contractions are cyclic, the 
driving force must satisfy the steady-motion constraint: < AF(i) > c = T c _1 J^ c AF(t)dt = 0, 
where T C (V) is the one-step ground-contact period (shown in Fig. 1). 

Assuming the displacements Ar to be small, we introduce the librational part of the 
potential energy AU[Ar(t)] in the harmonic approximation: AUq = ko(Ax 2 + Ay 2 )/2 
through the body stiffness coefficient fco(V). Combining this with the kinetic energy 
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AK = m(A x +A y )/2 and employing Lagrangian formalism, one deduces the New- 
ton equations mA r +A; Ar = 0. The free COM motion is therefore a superposition of the 
two linear oscillations: 

7T 7T 

Ax (t) = Al cos(uj t - -); Ay (t) = Ah sm(uj t - -). (1) 

These are solutions Ar (t) given by the harmonic amplitudes Al (V) and Ah (V) and the 
one-step angular frequency lj (V), with lu = 2n/T c = ^Jh^Jm. The backward elliptical 
COM rotation in Eq.(l) is due to the harmonic part of the driving force AFo(t) = — A;oAro, 
derived from the experiment and treated as an inertial force. Its components: 

AF 0x (t) = —mLulAl sm(ujot); AF 0y (t) = mcc^A/io cos(LU t), (2) 

are shown by solid lines in Fig. 1. With increasing speed, anharmonic displacements become 
important and therefore AU = AU + AU\. Without loss of generality, the anharmonic part 
of the mechanical potential energy AU\ is parametrized in terms of the anharmonic force 
amplitudes Al\ and A/ii, namely 
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AC/i(Ar) / ''° 



A/i 



Ah Ax 2 Ay + ^-Ay 3 + 0(A:rAy 2 ) + 0(Ax 3 ) 



(3) 



A/ 3A/i 

Within the perturbation scheme, the nonlinear forces are defined by the derivatives AF 1 = 
—dAUi/dAr taken at Ar = Ar given in Eq.(l). This results in 

F x (t) = —mul{Al sm(uj t) + Al 1 sm(2u t)}; 

F y (t) = mg + mujQ[Aho cos(uj t) — Ahi cos(2cj t)]- (4) 

The third and the fourth terms in Eq.(3) correspond to the terms of O[cos(2c;ot)] and of 
O[sin(2ujot)}, which formally should appear in F x and F y , respectively. Both the terms are 
omitted in Eqs.(4) because they are not observed in the available data on human walking. 
[7,8] Also, the steady-motion constraint provides the force-amplitude relation A/ A/i = 
A/i A/ii that is nevertheless violated, even in the case of the small V. As a matter of fact, 
the theory behind this relation presumes that AF 1 is a conservative force, which disagrees 
with the experimental data. 

Let us introduce a generalized velocity-dependent Lagrangian [9] AL(Ar,A r) = AK — 
AU e ff, where AC/ e // = AUo + AU\ — AK\ and AKi is the anharmonic kinetic energy. 
[10] Within the scope of this analysis, Eqs. (4) are not altered and the steady-motion 
constraint is satisfied by new kinetic terms. With the help of the frictional coefficient j(V), 
we also introduce the resistance force AF res (t) = —7 A r 1; associated with the anharmonic 
displacements Ar! = Ar — Ar . The latter obey the equations 

mA r\ +7 A r"i +A; Ari = AFi(t), (5) 

where AFi(t) is given by the last terms in Eqs. (4). Solutions of inhomogeneous differential 
equations (5) provide the desired description for the body's COM motion in the ground 
coordinate system: 

A/ x sin(2w ^ + V 9 ) 



X(t) = Vt + Ax (t) 



3 v/l+tan^) 2 ' 



Y{t) = H + Ay o{ t) + ^^tf±^. (6) 

o a/1 +tan(v?) i 
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They are found [9,10] for the steady motion regimes in the weak friction approximation, 
with 

<p(V) = arctan < -■ Ul (V) = 1. (7) 

As seen from Eqs.(6), these regimes are established not only by the frequencies, but also by 
the amplitudes. The latter statement follows from the force-time fitting analysis given in 
Fig. 1 that can be explicit in the force-amplitude ratio: A/o^/A^ 6 ^ = Ah^^ / Ah^^ 
which equals 2. As can be recognized from Eqs. (1), (6) and the force-amplitude ratio, a 
trajectory of the human body's COM, in the moving inertial coordinate system, is a closed 
orbit given by a shortened hypocycloid (Ari < Ar <C H), passing through three turning 
points in the backward direction. [11] For qualitative analysis, we describe this closed orbit 
by the characteristic ellipse, which crosses the same turning points and is introduced by the 



axes: Al = A/ + A/i/3a/1 + tan(y?) 2 and Ah = Ah + Ahi/3^/l + tan(y?) 2 , in the forward 
and vertical directions, respectively. One infers that the COM moves on the height H with 
the speed V, and simultaneously rotates along the hypocycloid circumscribed by a shrunken 
(or a flattened) ellipse of eccentricity e + (or e_), with e±(V) = a/1 — (Ali/Ahi) ±2 . 

The step-cycle external work W to t is performed by the COM to maintain the forward 
speed V and the height H, thus, W to t = Wf + W v . This work can be estimated through the 
instant power averaged over the cycle period: W to t(V) = 2nu <Wtot (t) > c , where Wtot 
(t) = F- R. Nevertheless, not all the components of the ground-reaction force, produced by 
active and passive muscles, contribute to the cyclic work. Bearing in mind the conditions i) of 
the periodicity of the driving force, ii) of the orthogonality between the linear and nonlinear 
displacements, and iii) of the conservative nature of the harmonic force, one deduces that 
the only nonzero contribution to W to t is due to the anharmonic part of the power Wi (t) = 
AFiA r 1 . This power follows from Eqs. (6) and provides 

Air ^ojAlj + Ahl) 
WtotiV) ~ T 1 + (2 7 /3mc ) 2 • (8) 

On the one hand, the total external work corresponds to that part of the mechanical energy 
that is lost as heat. It must be therefore restored in the next step through chemical energy 



by oxygen consumption. On the other hand, the external work is realized through the 
positive and negative contributions: W tot = W + — \W~\. These can be exemplified by the 
accelerated and decelerated forward body's displacements, respectively. A special case is 
ideal oscillational motion, when the two contributions are equal and, thus, W to t = 0. Also, 
the recovery coefficient, defined [3] for arbitrary cyclic motion as r = (W + — W to t)/W + , 
equals one. To estimate r(V), we specify the positive work by W + (V) = T C < P\W\ (t)\ >« 
where the auxiliary function P(x) = xH(x), with H(x) is the Heavyside step function. This 
leads to the recovery coefficient 

(p[cos(2u Q t) sm(2u t + <p)]) 

r(u , u x ) = 1 - ^ _ J . (9) 

2 v /(3^ /2cu 1 ) 2 + l 

If one employs the human- walk data on Uq XP \v) [12] and r^ exp \V), [3] Eq.(9) can be read 
as r(u^ xp \ uji) = r^ exp ^ and solved for the mass-specific frictional coefficient (see Fig. 2). 
Additionally, a straightforward estimation of the positive work performed in the forward 
(Wf~) and vertical (W+) directions, yields the anharmonic amplitudes 



Wt(l-r) 9o;2 + 4^2 / w + 

&h{V) = \j S + - ' ^(V)^^. (10) 

As follows from a numerical analysis of Eqs.(10) given in Fig. 3, the instability of the 
vertical COM librations occurs at the maximum speed for human walking with V max = 3.4 
m/s, in accord with the recent experimental data [5] vdax P ' ) = 3.2 m/s. The eccentricity of 
the orbit-characteristic ellipses (see Fig. 4) specifies this instability by the critical condition 
e -(Knax) — 1 associated with a dynamical transition from walking to running. There is also 
a dynamical crossover at the speed V cr = 1.7 m/s, which separates slow walking from fast 
walking. At this speed, e±(V cr ) = and the shrunken ellipses (A/ < Ah) transform into the 
flatter ellipses (Ah < Al). In reality, the crossover in locomotion is attributed to changes in 
performance of the human legs in the stiff-legged (Al < Ah) and the compliant (Ah < Al) 
walking. These two walking gaits are distinguished through distinct postures (and reaction- 
force records) and identified with a modern human walk and a walk of nonhuman primates, 
respectively. [5] 
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In conclusion, we have discussed the problem of animal locomotion on level ground 
in view of human data on reaction force and external work for walking at steady speeds. 
Application of the standard Lagrangian formalism developed for nonconservative forces, 
permits one to introduce a generalized equation for animal locomotion in the point-body- 
mass approximation. Such a description involves only integrative properties of body, given 
through the inertial, elastic, and resistance characteristics, and is therefore expected to 
be helpful for comparative studies of quadruped trotting of animals where lateral effects 
are also negligible. For the case of human walking, the locomotion for the body's COM 
is a superposition of the rectilinear motion with its backward rotation along a shortened 
hypocycloid. 

Records of the ground-reaction force elucidate a variety of body's functions. In animal 
locomotion, the observed reaction force acts as a motor-brake force, which additionally 
supports the body weight and controls the stability of forward advancement. The efficiency 
of the employed above body-support function is restricted by anatomical adaptation of 
the long-bone limbs for peak force-ground contacts. This adaptation was proven [13] to 
be universal for all terrestrial mammals, with evolution of their body mass. As to motor- 
brake effects, we have deduced that the velocity-independent forces, which produce harmonic 
body vibrations, are responsible for an effective-exchange mechanism of elastic mechanical 
energy, attributed to passive muscles and parametrized by the body stiffness coefficient. 
The nonlinear anharmonic librations, which are due to the nonconservative part of the 
driving force, produce the main part of the external work. The feedback between linear and 
nonlinear COM librations is revealed through the phase, the amplitude and the frequency 
constraints imposed on cyclic human walking. Finally, our analysis of the external work 
given without recourse to a swing-pendulum or spring-mass modeling provides insights into 
the two principal gaits of biped walking. We have seen that the stiff-to-compliant crossover 
in human walking arises from the changes of body-resistance performance and not from the 
body-elastic adjustment, prescribed earlier by the pendulum dynamics. 
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FIGURE CAPTURES 

Fig. 1. Analysis of force-platform data on human level walking. Records for the hori- 
zontal (F x ) and vertical (F y ) forces are taken from Fig. 1 in Ref. [8] (the case of mass 59 
kg and of speed 3.9 km/h) and given for the two-step period 2T C . Solid lines describe the 
force-time linear harmonics at distinct phases (denoted by letters) explicit in Eq.(2) along 
with the fitting amplitudes AZ = 0.012 m and Ah = 0.016 m. The lines indicated by 
hatched areas correspond to the double-frequency harmonics given in Eqs.(4) and adjusted 
with AZi = 0.006 m and Ah x = 0.008 m; H (« Im) stands for the body's COM height. 
Asymmetric deviations are due to the differences between i) the exertions by left and right 
feet and ii) the time intervals for the single-foot (T sc ) and double-foot (T dc ) ground contacts, 
with T c = T sc + T dc . Rotations of the driving force and the corresponding velocity are shown 
by dotted lines in the coordinates AF y (t) vs AF x (t) and AV y (t) vs AV x (t), respectively. 

Fig. 2. Human- walk characteristic frequencies against steady speed. Points for the 
cyclic frequency uq{V) reproduce cinematographic data (open circles) taken from Fig. 3 in 
Ref. [12]. They are extended by the case analyzed in Fig. 1 (open square) and fitted by 
Vo(V) = 4.94 + 4.02V (solid line). Points for the mass-specific frictional coefficient r y(V)/m 
(shaded circles) are found through Eq.(9) with the help of ujo(V) and r^ exp ^(V) (reproduced 
in the insert from Fig. 2 in Ref. [3]) and fitted by wi(V) = 6.37 - 6.15V + 2.38V 2 (given by 
the solid curve). 

Fig. 3. Anharmonic amplitudes for the driving force against speed in human walking. 
The horizontal and vertical amplitudes Ali(V) and Ahi(V), reduced by AZ (1.1) = 0.012 
to, are given by open and shaded circles, respectively. They are estimated through Eqs.(10) 
with the help of the parameters obtained in Fig. 2 and of data on the positive work W^ cxp ^ 
and Wv^^ (taken from Fig. 2 in Ref. [3]). The positive forward and vertical works are 
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adjusted at V = 1.1 m/s (shown by open and shaded squares, respectively). The curves are 
third-order polynomial fits extrapolated to the maximum speed V max = 3.4 m/s (shown by 
the arrow). 

Fig. 4. The human center-of-mass body orbit characteristics for slow and fast level 
walking. Inserts show ellipses, which circumscribe the closed COM trajectories in the inertial 
coordinate system moving with different speeds V (indicated by arrows). The elliptic axes 
are reduced by the amplitude Al = 0.012 m. Experimental data on the elliptic eccentricities 
(shown by open circles) for the shrunken and flatter ellipses obtained, respectively, through 
yjl - W+ (cxp) /W v +(cxp) and ^/l - W v Hcxp) /Wp cxp) , with the help of the experimental data 
reported in Ref. [3]. The fitting curve corresponds to the analysis given in Fig. 3. 
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